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PHASE TRANSITIONS FOR MODIFIED ERDOS RENYI 

PROCESSES 

SVANTE JANSON AND JOEL SPENCER 

(N 

^V Abstract. A fundamental and very well studied region of the Erdos- 

Renyi process is the phase transition at m ~ ^ edges in which a giant 
component suddenly appears. We examine the process beginning with 
an initial graph. We further examine the Bohman-Frieze process in 
which edges between isolated vertices are more likely. While the posi- 
tions of the phase transitions vary, the three processes belong, roughly 
speaking, to the same universality class. In particular, the growth of the 
(**N , giant component in the barely supercritical region is linear in all cases. 

U' 

1. Introduction 

( n ) 
The standard Erdos-Renyi process (G(n,m))^L starts with an empty 

graph G(n, 0) = E n with n vertices and adds edges one by one in ran- 

£> ■ dom order, uniformly over all possibilities, i.e., drawing the edges uniformly 

without replacement. (Hence, G(n,m) has n vertices and m edges.) This 

^Li_ | random graph model has been studied a great deal, starting with Erdos and 

t^J- ■ Renyi [g, lZ(, see for example the monographs by Bollobas [4j] and Janson, 

Luczak and Rucihski [12]. 

O ' The purpose of this paper is to study two modifications of this process. We 

are interested in the sizes (orders) of the components of the random graphs; 
in particular whether there exists a giant component of size comparable to 
the entire graph and, if so, how large it is. (We ignore the internal structure 

K^ , of the components.) We denote the components of a graph G by Ci(G), 

i = 1, . . . ,v(G), where thus v(G) is the number of components of G, and 
their sizes by Q(G) := |Cj(G7)|, 1 < i < v(G); we will always assume that 

the components are ordered such that C\ > C2 > (For convenience we 

also define Ci{G) = when i > v(G).) We will often, as just done, omit the 
argument G when the graph is clear from the context. We further denote 
the edge set of G by E(G), the number of edges by e(G7) := |_E(G7)|, and the 
number of vertices by |67| (the order or size of 67). 

We recall the fundamental result for G(n, m) [7( that if n — > 00 and m ~ 
cn/2 for some constant c, then C\ = p(c)n + o p (n), where p(c) = if c < 1, 
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2 SVANTE JANSON AND JOEL SPENCER 

and p(c) > if c > 1. (Furthermore, C2 = o p (n) for every c.) This is 
usually expressed by saying that there is a threshold or phase transition at 
m = n/2. See further @; 0; [H]. Moreover, as <5 \ 0, p(l + (5) ~ 25 (see [1, 
Theorem 3.17] for a generalization to certain other random graphs). (For 
the notation OptnV and other standard notations used below such as w.h.p., 
see e.g. [H] and Q.) 

In the first modification of the Erdos-Renyi process, we assume that some 
(non-random) edges are present initially; additional edges then are added 
randomly as above. We actually consider three slightly different versions of 
this process; see Section [2] for details. Our main result for these processes 
(Theorem I2.ip characterizes the existence and size of a giant component 
in terms of the initial edges (more precisely, the sizes of the components 
defined by them) and the number of added random edges. We define the 
susceptibility S2 as the average size of the component containing a random 
vertex in the initial graph, see ()2.ip - (j2.3p . and show the existence of a 
threshold when t c n/2 edges are added, where t c := s^ . (This was also done, 
under a technical assumption, in Spencer and Wormald |l4|].) Moreover, we 
give upper and lower bounds for the size of the giant component after the 
threshold in terms of S2 and two related quantities (higher moments of the 
component size) S3 and S4 for the initial graph, also defined in (|2.ip - (|2.3p . 

Our second modification is known as the Bohman-Frieze process, after 
Bohman and Frieze [2J. The initial graph on n vertices is empty. At each 
round two edges e% = {vi,w\} and e2 = {^2,^2} are selected indepen- 
dently and uniformly. If both v\ and w\ are isolated vertices the edge e\ 
is added to the graph; otherwise the edge e2 is added to the graph. We 
let BF m denote this process when m edges are added. This is a natural 
example of an Achlioptas process, in which a choice may be made from two 
randomly chosen potential edges. In Bohman and Frieze |2j and Bohman, 
Frieze and Wormald [3j it was shown that the phase transition is deferred 



beyond m ~ n/2. More precisely, it is proved in Spencer and Wormald 14] 
that the Bohman-Frieze process has a phase transition at some t c « 1.1763. 
In the present paper we study further what happens just after the phase 
transition, using the result just described for the Erdos-Renyi process with 
initial edges. The idea is, as in [lj], that to study the process at a time 
£1 > t c , we stop the process at a suitable time to just before the phase 
transition, and then approximate the evolution between to and t\ by an 
Erdos-Renyi process, using the graph obtained at time to as our initial 
graph. In order to apply Theorem 12.11 we then need information on s?, S3 
and S4 in the subcritical phase. The analysis in Spencer and Wormald pjj] 
of the Bohman-Frieze process (and a class of generalizations of it) is based 
on studying the susceptibility S2 in the subcritical region. We will use some 
results from [14J] , reviewed in Section [3l and extend them to S3 and S4 in 
order to obtain the required results needed to apply Theorem 12.11 

In particular, we show that after the phase transition, the giant compo- 
nent grows at a linear rate, just as for the Erdos-Renyi process. The precise 
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statement is given by Theorem 13.51 The original Erdos-Renyi process, the 
process from an appropriate starting point, and the Bohman-Frieze process 
appear to be in what mathematical physicists loosely call the same univer- 
sality class. While the placement of the phase transitions differ the nature 
of the phase transitions appears to be basically the same. A very different 
picture was given for a related process in [1]. There, as in the Bohman- 
Frieze process, two random potential edges e\ = {v\,w{\ and &2 = {^2,^2} 
are given. However the edge is selected by the Product Rule: we select that 
edge for which the product of the component sizes of the two vertices is 
largest. Strong computational evidence is presented indicating clearly that 
this process is not in the same univerality class as the three processes we 
compare. We feel, nonetheless, that there is likely to be a wide variety of 
processes in the same universality class as the bedrock Erdos-Renyi process. 

The main results are stated in Sections [2] and El and proved in Sections 0] 
and[5j 

Our results are asymptotic, as the size grows. All unspecified limits are 
as n — > 00. 

Acknowledgement. This research was mainly done at Institute Mittag- 
Leffler, Djursholm, Sweden, during the program Discrete Probability, 2009. 
We thank other participants, in particular Oliver Riordan, for helpful com- 
ments. 

We thank Will Perkins for the numerical calculations in Remark 13.61 

2. ERDOS-RENYI PROCESS WITH AN INITIAL GRAPH 

The purpose of this section is to study the Erdos-Renyi process when 
some edges are present initially. We define three different but closely related 
versions of the process. 

Let F be a subgraph of K n with vertex set V(F) = V(K n ) = {1, . . . , n}. 

Define (G(m, n;F))^i. by starting with G(n,0;F) := F and adding 
the (2) — e(F) edges in E(K n ) \ E{F) one by one in random order, i.e., by 
drawing without replacement. 

For our purposes it will be convenient to consider two modifications of this 
random graph process. (Both modifications are well-known for G(n,m).) 
We define (G*(n, m; F))™ =Q by starting with G*(n,0;F) := F and then 
adding at each time step an edge randomly drawn (with replacement) from 
E(K n ), provided this edge is not already present (in which case nothing 
happens). In particular, G*(n, m) := G*(n, m; E n ) is defined as G(n, m) but 
drawing the edges with replacement. In general, we have E(G*(n,m;F)) = 
E(G*(n,m))UE(F). 

Note that the number of edges in G*(n,m) may be less than m. Al- 
ternatively, we may regard G*(n,m;F) as a multigraph and add the edges 
whether they already are present or not; then the number of edges is always 
exactly m + e(F). Since we will study the component sizes only, this makes 
no difference for the present paper. 
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The second modification is to use continuous time. We may think of the 
(2) edges as arriving according to independent Poisson processes with rates 
1/n; thus edges appear at a total rate Q) fn = ^p and each edge is chosen 
uniformly at random and independently of all previous choices. We define 
G(n,t;F) to be F together with all edges that have arrived in [0, t]. (As 
above, we can consider either a multigraph version or the corresponding 
process of simple graphs, obtained by ignoring all edges that already appear 
in the graph.) Hence, if i and j are two vertices that are not already joined 
by an edge in F, then the probability that they are joined in G(n, t; F) is 
1 — e~ l ' n = t/n + 0(t 2 /n 2 ), and these events are independent for different 
pairs i,j. (Starting with the empty graph we thus obtain G(n,p) with p = 
\—e~ t ' n . We could change the time scale slightly to obtain exactly G(n, t/n), 
and asymptotically we obtain the same results for the two versions.) 

Note that if N(t) is the total number of edges arriving in [0,£], then 
N(t) ~ PofQ) t/n) = Po(^2— t), and, with an obvious coupling of the pro- 
cesses, G(n,t;F) = G*(n, N(t); F). For constant t, N{t)/{n/2) -£-> t as 
n — > 00 by the law of large numbers. Moreover, the expected number of 
repeated edges in G*(n,m;F) is at most (^/(J) +™>\E(F)\/(£); if for ex- 
ample, as in Theorem 12.11 m = 0(n) and |i?(.F)| = O(n), then this is 0(1), 
which will be negligible. Standard arguments, comparing the processes at 
times t and (l±e)t, show that for the properties considered here, and asymp- 
totically as n — > 00, we then obtain the same results for G(n, \nt/2\;F), 
G*(n, [nt/2\;F), and G(n,t;F). 

We define, for a graph G with components of sizes C±, . . . , C v , and k > 1, 

S k = S k (G):=Y,Ct (2.1) 



summing over all components of G. Thus S\(G) = \G\, the number of 
vertices. We normalize these sums by dividing by \G\ and define 

Hence, s\(G) = 1 for every G. Note that 

Sk(G) = Y,§- ] Ct\ (2.3) 

i 

which is the (k — l):th moment of the size of the component containing 
a randomly chosen vertex. In particular, S2(G) is the average size of the 
component containing a random vertex. The number 52(G) is called the 
susceptibility; see e.g. [Ill; ll3j; llO| for results on the susceptibility in G(n,m) 



and some other random graphs. 
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It follows from the definitions (|2.ip and (|2.2[) that Sk and Sk are (weakly) 
increasing in k; in particular, Sk(G) > si(G) = 1 for every A; and G. More- 
over, Holder's inequality and ()2.3[) imply that the stronger result that s fc 

(and even s fc , k > 2) is (weakly) increasing in fe. 

Note further that the number of edges in a component of size C- h is at 
most ( 2 ') < C?; hence, for any graph G, 

\E(G)\ < S 2 (G). (2.4) 

We will use these functionals for the initial graph F to characterize the 
existence and size of a giant component in the random graph processes 
starting with F. An informal summary of the following theorem (our main 
result in this section) is that there is a phase transition at t c := 1/ 82(F), 
and that for t = t c + 5 with 5 small, there is a giant component of size 
« 2(s2{F) 3 j ss(F))Sn. For the special case when F = E n is empty, s 2 = 
S3 = 1 and we recover the well-known result for the Erdds-Renyi process 
mentioned above that there is a phase transition at t c = 1 (i.e., at re/2 edges) 
and further for t = 1 + S, there is a giant component of size « 2<5re. The 
formal statement is asymptotic, and we thus consider a sequence F n . 

Theorem 2.1. Suppose that for each n (at least in some subsequence), F n 
is a given graph with n vertices, and suppose that sup n S3(F n ) < 00. Let 
the random variable Z n be the size of the component containing a random 
vertex in F n . 

Consider the random graph processes G(n,t;F n ). Then, for any fixed 
t > 0, the following hold as n — )• 00, with st '■= Sk{F n ), 

(i) Ift < l/8 it then Ci(G(n,t;F n )) = o p {n). 
(ii) If t > l/s2, then there is a unique p n > such that 

p n = l-Ee" pA , 

and we have 

C\(G{n,t;F n )) = p n n + o v (n). 

(hi) Ift > l/s 2 , let 5 n :=t- l/s 2 > 0. Then 

re s 3 

If further <5 n s i s 4/ s i — §> then also 

d(G(re,i;F n )) ^ ^ s| ^ , 8^ s|s 4 

4 



<25 n ^[l + -5 n ^\+o p (l). 
n s 3 V 3 si ' 



(iv) In (hi), if in addition lim inf n _>oo ^n > 0, then moreover w.h.p. 

C^Gjn^Fn)) 4 

> 2d„— (1 - 2d n S2) 

n s 3 
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and, if <5 n s|s 4 /s| < §, 



Ci(G(n,t;F n )) /W ^/8 J s|s 4 



< 1&n- 1 + ^n — 

n S3 \ 3 S3 

TTie same results hold for the random graph processes G(n, [nt/2\ ; F) and 
G*(n,[nt/2\;F). 

The proof is given in Section [H Note that by (|2.3|) , 

EZ* = s k +i(F n ), k>\. (2.5) 

3. The Bohman-Frieze process 

Recall the definition of the Bohman-Frieze process from Section [IJ see 
[2j; y; ll4|: we are at each round presented with two random edges e\ = 
{vi, W\} and ei = {t>2, 1^2} in the complete graph K n and choose one of them; 
we choose e\ if both its endpoints v% and w\ are isolated, and otherwise we 
choose e2- We let BF m denote the random graph created by this process 
when m edges are added. (The size n is not shown explicitly.) We further 
define, using the natural time scale, BF(t) := BFi nt / 2 j . (For convenience, we 
sometimes omit rounding to integers in expressions below.) 

Note that if we add e%, then it always joins two previously isolated vertices, 
while if we add &i, it is uniformly distributed and independent of the existing 
graph. We call the added edges ei Erdos-Renyi edges, since all edges in the 
Erdos-Renyi process are of this type. 

Remark 3.1. We have talked about edges e± and e^, but it is technically 
convenient in the proofs to allow also loops (as in [14| ) ; we thus assume in the 
proofs below that in each round, the vertices V\,iU\,Vi-,Wi are independent, 
uniformly distributed, random vertices. It is easily seen that the results 
proved for this version hold also if we assume that there are no loops, for 
example by conditioning on the event that no loops are presented during 
the first nt/2 rounds; we omit the details. 

For a graph G, let rii = rii(G) be the number of vertices in components of 
order i, and let Xi = Xi(G) := m(G)/\G\ be the proportion of the total num- 
ber of vertices that are in such components. (Thus, Sk{G) = ^2ii k ~ 1 Xi(G).) 
For the Bohman-Frieze process, we need only m, the number of isolated 
vertices, and the corresponding proportion x\ := n\jn. 

For the Bohman-Frieze proce ss ( and some generalizations of it), it is 
shown in Spencer and Wormald [lJ] that the random variables xi(BF(t)) 
(for any fixed t < 00) and S2(BF(t)) (for any fixed t < t c ) converge in 
probability, as n — >• 00, to some deterministic values x\{t) and s~2(£); these 
limit values are given as solutions to differential equations. We extend this 
to S3 and S4 as follows. 



We first define, as in [14|], the deterministic function x\{t) as the solution 
to the differential equation 

x ' 1 (t) = -xl(t)-(l-x 2 l (t))xi(t), t>0, (3.1) 
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with initial condition xi(0) = 1; by [14! . Theorem 2.1], x\(t) is defined and 
positive for all t > 0, and by [14J, Theorem 1.1], xi(BF(t)) — > x\(t) for 
every fixed t > 0. 

We further define functions s~2(t), s 3 (t), s^t) as the solutions to the dif- 
ferential equations 

s' 2 (t) = x 2 1 (t) + {l-x 2 1 (t))s 2 2 (t), (3.2) 

s' 3 (t) = 3xl(t) + 3(1 - xl(t))s 2 (t)s 3 (t), (3.3) 

s' A (t) = 7x\{t) + (1 - x\{t)) (4s 2 (t)s 4 (t) + 3fl(t)) , (3.4) 
with initial conditions 

52(0) = s 3 (0) = s 4 (0) = 1. (3.5) 



The function s~2(t) is studied in Spencer and Wormald Jjj, Theorem 2.2], 
and it is shown there that it explodes at some finite t c , i.e., the solution s~2{t) 
is (uniquely) defined for t G [0, t c ), but s~2(t) /* +oo as t / i c ; it is further 
shown 1J, Theorem 1.1] that this t c is the time of the phase transition for 
the Bohman-Frieze process, when a giant component first appears, and that 
for any fixed t < t c , s~2(BF(i)) — > S2(t). We extend these results to S3 and 
S4 as follows. 

Theorem 3.2. The functions s~2(t), s 3 (t), s^it) are uniquely defined by 
(|3.2p - (|3.5p for all t G [0, t c ). As t f- t c , there exist positive constants a and 
/3 such that 

, . a 

s 2 (t) 



t c -t 
s 3 (t) ~ PUtf 



/3a 3 



2^ U\5 



s 4 (t) ~ 3(3 2 s 2 (t) 



(tc - i) 3 ' 

3/3 2 a 5 



More precisely, s k (t) = a k (t c - i)-( 2fe ~ 3 )(l + 0(t c - t)) for k = 2,3,4 with 
a 2 = a, 03 = /3a 3 , 04 = 3/3 2 a 5 . 

W^e /taue a = (l — x 2 (£ c )) , w//«te /3 = <?(£ c ) *s gwen by (|5.2p and (|5.5p . 

Theorem 3.3. For any fixed t £ [0,t c ), and k = 2,3,4, Sfc(BF(i)) — >Sk{t). 

Remark 3.4. It is straightforward to extend Theorem 13.31 to any k > 2, 
with s~k(t) given by a differential equation similar to (|3.2|) - (|3.4|) (involving 
Sj for j < k, so the functions are defined recursively). We leave the details 
to the reader since we only use k < 4 in the present paper. 

Proofs are given in Section [5] Using these results for the subcritical phase, 
we obtain the following for the supercritical phase; again the proof is given 
in Section [5j 
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Theorem 3.5. There exists constants 7 = 2(1 — x\ (t c ))/(3 > and K < 00 
such that for any fixed 5 > 0, w.h.p. 



KlV >^m*.+ii)) KtV ,_ 



7<5 

n 

Remark 3.6. Numerical calculations of Will Perkins give t c ~ 1.1763, 
xi(t c ) « 0.2438, a » 1.063, /? « 0.764, a 2 = a, a 3 ^ 0.917, a 4 » 2.375 
and 7 « 2.463. 

There is an obvious conjecture (made explicit in [14]]) that Ci(BF(t))/n — )• 
PBf(^) for some function pbf : [0, 00) — > [0,1]; equivalently, Ci(BF(£)) = 
PBf(t)n + o p (n). (For £ < t c , clearly this holds with PBp(t) = 0.) In Spencer 
and Wormald [14] it was further conjectured that lim^ t + Pbf(^) = 0; in the 
language of Mathematical Physics, this says that the phase transition is not 
first order. If such an /9bf exists, Theorem 13.51 resolves the latter conjecture 
positively and further gives the asymptotic behavior pBF(t c + $) ~ 7^ as 
5^0+. 

Remark 3.7. We further conjecture that the function p&f is smooth on 
[t c ,oo); if this is the case, then Theorem 13.51 shows that p BF (i^) = 7. This 
conjecture would imply that 5 4 ' 3 in Theorem 13.51 could be replaced by S 2 ; 
unfortunately, our approximations are not sharp enough to show this. 

4. Proof of Theorem 12.11 

We begin with a simple lemma (related to results in [5J, Section 5]). 

Lemma 4.1. Let Y > be a random variable with 1 < E7 < 00. 
(i) There is a unique p > such that 

p=l-Ee- pY . (4.1) 

(ii) If EY 2 < 00, then 

2(EY-1) 

P> EY 2 ■ 
(iii) If EY 3 < 00 and 8(EY -1)EY 3 < 3(E7 2 ) 2 , then 



3EY 2 - a/9(E Y 2 ) 2 - 24(E Y - 1)EY 3 
P< 2EY 3 

4(E7-1) 



EY 2 + J(EY 2 ) 2 - |(Ey - 1)E7 3 
2(EY-1) / 8(Ey-l)Ey 3 



EY 2 V 3(Er 2 ) 2 

(iv) Let Y n , n > 1, be random variables with Y n > and EY n > 1 and 
let p n > be the corresponding numbers such that p n = 1 — Ee~ pnYn . If 

Y n — > Y for some Y with E7> 1, then p n — > p > satisfying (|4.ip . On 

the other hand, if Y n — > Y with EY"< 1, then p n — > 0. 
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Remark 4.2. In fact, (14. ip is the standard equation for the survival prob- 
ability of a Galton-Watson process with a mixed Poisson Po(Y) offspring 
distribution. Parts (i) and (iv) follow easily from standard results on branch- 
ing processes. We prefer, however, to give direct proofs (also easy). Note 
further that p = always is another solution to (|4.ip . If EY < 1, then 
p = is the only non-negative solution, either by branching process theory, 
or because 

1 - E e~ sY = E(l - e~ sY ) < E(sY) < s (4.2) 

for every s > 0, with strict inequality unless sY = a.e. and E(sY) = s, 
which together imply s = 0. 

Proof. The function <p(s) := 1 — E e~ , s E [0, oo), is increasing and concave 
with < <p(s) < 1, <p(0) = and <p'(0) = EY > 1. Consequently, <p(s) > s 
for small s > 0, but <p(s) < s for s > 1, say, and there is a unique p > 



such that </?(p) = p. This proves (i) Note that (p(s) > s for < s < p and 
tp(s) < s for s > p. 



We next prove (iv) If E Y > 1, let < e < p. Then </?(p — e) > p — e and 



thus, because Y n — > Y, 

1 - E e -^- £ ) y ™ ->. 1 - E e -^ £ ) y > p - e, 

so for large n, 1 — Ee _ ' /3_e - ) n > p — e and thus p — e < p n . Similarly, for 
large n, 1 — Ee - '^' 5 " < p + e and thus p + e > p n . Since e is arbitrarily 
small, it follows that p n —> p. 

If instead EY < 1, then <p( s ) < s f° r every s > by (|4.2p and the 
comment after it. Hence the same argument shows that for every e > 0, 
p n < e for large n; thus p n — >• 0. 



To see (ii) , observe that e x < 1 — x+x 2 /2 for x > 0, with strict inequality 
unless x = 0, and thus, when EY 2 < oo, 

/ n 2 Y 2 \ n 2 

p = E(l-e- pY ) >E(pY-^ J =pEY-^-EY 2 . 



Hence, 1 > E Y - pE Y 2 /2, which yields (ii 



For (hi) , we first note that, similarly, e x > 1 — x + x 2 /2 — x 3 /6 for x > 0, 
again with strict inequality unless x = 0, and thus, provided EY 3 < oo, 

/ n 2y2 n 3v 3 \ n 2 n 3 

p = E(l-e-P Y ) <E(pY-^^ + ^^) =/0 EY-^EY 2 + ^EY 3 . 
V 2 6 / 2 6 



This can be written 

EY 3 p 2 -3EY 2 p + 6(EY-l) > 0. (4.3) 

As long as the discriminant 9(EY 2 ) 2 - 24(EY - 1)EY 3 > 0, the corre- 
sponding quadratic equation (with equality instead of >) has two roots 

_ 3EY 2 ± ^9(1 Y 2 ) 2 - 24(E Y - 1) E Y 3 
P± ~ 2EY 3 
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and we have either p < p_ or p > p + . In order to rule out the latter 
possibility, we consider the random variable Yt := tY for to < t < 1, where 
to = 1/EY. Note that for t < t < 1, E Y t > 1 and thus there is an p(t) > 
such that p(t) = 1 — Ee - ' 1 "' '; by (iv) , p(t) is a continuous function of t. 



Further, for to < t < 1, 

9(EY ( 2 ) 2 = 9t 4 (EY 2 ) 2 > 24t 4 (Ey-l)Ey 3 = 24(Eyi-t)Ey t 3 

> 24(EY t - l)Elf; 

hence the discriminant is non- negative for each Yt, and there are corre- 
sponding roots p±(t). These are continuous functions of t and for each 

t G (t ,l), p(t) < p-(t) or p{t) > p + (t). As t \ t , El* -» 1 and 

p+(t) -> 3E^/(2El f 3) > o while, by [(iv)] again, p(t) -»• 0. Hence, 



pit) < p+(t) for t close to to, and by continuity, p(t) < p+(t) for all t G (to, 1] 
(since equality is impossible by (|4.3j) ). Consequently, p < p + and thus 
p< p-. 

Finally, we use straightforward algebra and the fact that for x £ [0,1], 
\/l — x > (1 — x)/(l + x) and thus 

1 1 + x 

< — =— . □ 



i + yr^x - 2 

Proof of Theorem \2.1\ Note that the assumptions and (|2.4p imply that 

1^(^)1 < S 3 (F n ) = ns 3 (F n ) = 0{n). 

Hence, by the discussion in Section [21 it suffices to consider G(n,t;F n ). 

The main idea is that we may collapse each component Ci(F n ) of F n to a 
"supervertex" with weight 

x i = x\ n) :=\C i (F n )\=C i (F n ). (4.4) 

The probability of an edge between Ci(F n ) and Cj(F n ) in G(n,t;F n ) is, for 

Mil 

Pii (t) = l- e - te ^/ n . (4.5) 

Hence, to obtain the distribution of component sizes in G(n, t; F n ) we may in- 
stead consider the random graph H n with v = v{F n ) vertices having weights 
Xi given by (14. 4p and edges added independently with probabilities pij given 
by (I4.5J) : note that the size of a component in G(n,t;F n ) is given by the 
weight of the corresponding component in H n , i.e., the sum of the weights 
of the vertices in it. 

The random graph H n is an instance of the general random graph model 
studied in Bollobas, Janson and Riordan [5J]; we will use results from [5J, 
and therefore we show the relation in some detail. 

We will actually consider a subsequence only, for technical reasons, and 
thus we at first obtain the result for this subsequence only. However, this 
means that if we start with any subsequence of the original sequence, there 
exists a subsubsequence where the result holds; this fact implies that the 
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result actually holds for the full sequence by the subsubsequence principle, 
see e.g. [12J, p. 12]. 

We have defined Z n as the size of the component containing a random 
vertex in F n . Let v n be the distribution of Z n ; thus v n is the probability 
measure on Z+ := {1, 2, . . . } given by £\ %8 Ci ■ By ^B) , E Z n = s 2 (F n ) < 
s 3(Fn) = 0(1), which implies that the sequence of random variables Z n is 
tight, see e.g. [8, Section 5.8.3]. Consequently (see [8, Theorem 5.8.5]), we 
may select a subsequence such that Z n converges in distribution to some 
random variable Z. Equivalently, v n converges (weakly) to some probability 
measure ju on Z+, where u, is the distribution of Z. Moreover, EZ^ = 
s 3(Fn) = 0(1), and thus [8j, Theorem 5.4.2] Z n are uniformly integrable; 
consequently [8j, Theorem 5.5.8], S2(F n ) = ~E Z n — > EZ. We denote this 
limit by s 2 , and have thus 

s 2 (F n )^s 2 = EZ. (4.6) 

Let Vk(F n ) be the number of components of order k in F n and let u n be 
the measure on Z + defined by 

Mk} ■= ^ 

n 

Equivalently, u n := ^X^iLi^Ci- The total mass of u n is thus f n (Z + ) = 
v(F n )/n < 1. (In general, ^ n is not a probability measure.) 

The total size of the components of order k in F n is kvk(F n ), and thus 

£„{&} = F(Z n = k) = kVk{Fn) = ku n {k}. 

n 

Let \jl be the measure on Z + given by 

fi{k} := Jl{k}/k, k > 1. 

Since we have Pn{&} — >■ ju{&}, we also have 

v n {k] = D n {k}/k ->■ £{&}/& = //{&} 

for every A; > 1. Moreover, if / : Z+ — )• M is any bounded function, and 
p(fc) := f(k)/k, then the convergence z? n — )• /2 implies 

f(x)dv n (x)= I g(x)du n (x) ->• / g(x)dfi(x) = f(x)d/u,(x). 

Jz+ Jz+ JZ+ 

Hence f n — >• it weakly; in particular 

i/„(.A) -)■ //(A) for every ,4 C Z+. (4.7) 

We let (x n ) n >x be the sequence (Ci(F„), . . . , C„ n (F n )) of component sizes 
of F n , where v n := v(F n ). We have just shown that the triple V := 
(Z + , /i, (x ra )„>i) is a generalized vertex space in the sense of [5j, p. 10]; in 
particular, the crucial condition [5j, (2.4)] is our (|4.7p . 

We define the kernel k on Z + by 

k(x,h) := to/ (4-8) 
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(recall that t is fixed); the probability (J4.5)) of an edge in H n between (su- 
per)vertices with weights x% and Xj is thus 1 — exp(— K,(xi,Xj)/n), which 
agrees with ja, (2.6)]. Hence, our random graph H n is the graph denoted 
G v (u,k) in ||. 

We further have, with Xi = Ci(F n ), by (|4.5p . 



~Ee(iT n ) = - y~] Pij = - J~] (l - exp(-tXiXj/n)) 

1 v- i/l^ V t 



and 



„ oo oo 

/ xd/j 1 (x) = y^xdfi{x} = \^ djl{x} = ju(Z + ) = 1 
•^ z + z=l x=l 



(4.9) 



(since ju is a probability measure on Z + ); hence 

k(x, y) d/i(x) d/i(y) = t ( / xd^(x) I -t (4.10) 



and 

-Ee(H n )<- K(x J y)dfi(x)dfJ,(y). 

Together with [5|, Lemma 8.1], this shows that 

- E e(H n ) ->■ - / / K(x,y)d/j,(x)diJ,(y), 

and thus, using also (|4.10p . the kernel k is graphical [a, Definition 2.7]. 

We can now apply the results of [a]. The kernel K(x,y) is of the special 
type ip(x)ip(y) (with i/j(x) := i 1 ' 2 ^), which is the ranfc 1 case studied in [a, 
Section 16.4], and it follows by [5J, Theorem 3.1 and (16.8)] that H n has a 
giant component if and only if ||T K || > 1, where T K is the integral operator 
with kernel k; in the rank 1 case T K has the norm, using also (j4.6p . 



II^kII = / ip{x) djx{x) = j tx d/J,(x) = / txdjl(x) = t~EZ = ts2- 

Jz + Ji+ Jz + 

Hence there is a phase transition at t c := l/s2- We consider the cases 
t < t c and t > t c separately. 

4.1. The (sub)critical case. Consider first the case t < s^~ ; then H n thus 
has no giant component; more precisely, 

Ci(fln) = op(n). (4.11) 

Recall, however, that we really are interested in the size of the largest com- 
ponent of G(n, t; F n ), which is the same as the largest weight of a component 
in H n . (Note also that the component with largest weight not necessarily 
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is the component with largest number of vertices.) Nevertheless, the corre- 
sponding estimate follows easily: Let A > 0. Then the total weight of all 
vertices in H n of weight larger than A is 

J>il[xi > A] = Y, kv k (F n ) < A- x J2k 2 v k (F n ) = A~ l S 2 (F n ) 

i k>A k>l 

= nA- 1 s 2 (F n ), 

and thus the weight of any component C in H n is 

^Xi <^2xil[xi < A} + J2xil[xi > A] < A\C\ + nA~ 1 s 2 (F n ) 
iec iec i 

<AC 1 (H n ) + nA- 1 s 2 (F n ). 

For any e > 0, we may choose A = A n := e~ 1 s 2 (F n ) and find (since A n = 
0(1)) w.h.p., using (J37ITJ ). 

Ci(G(n,t;F n )) = sup^Xi < A n d(H n ) + en < 2en. (4.12) 

C i£C 

which proves (i) when t < l/s 2 . 

4.2. The supercritical case. Suppose now that t > s^ ■ 

By |5|, Theorem 3.1], the size Ci{H n ) of the largest component C\ of H n 
satisfies 

Ci{H n ) P 

> p{K) > 0. 

n 

Furthermore C 2 {H n ) = o p (n), and it follows by the same argument as for 

(|4.12p above that the weigth of any component C ^ C\ of H n is at most 

max) Xi < A n C 2 (H n ) + en < 2en 

^ iec 

w.h.p., and thus o p (n). Since C\ has weight > \C\\ = p(n)n + o p (n), it 
follows that w.h.p. the largest component C\ of H n also has the largest 
weight, and thus corresponds to the largest component in G(n,t;F n ), while 
C 2 (G(n,t;F n ))=o p (n). 

It remains to find the weight of C\. We first note that by [5J, (2.13), 
Theorem 6.2 and (5.3)], p(n) = L p K (x)dp(x), where p K {x) is the unique 
positive solution to 

p K = $„(/>«) :=l-e- T « p \ 
Since 

T K p K (x) := / K(x,y)p K (y)dp(y) = tx / yp K (y)dp(y), 



we thus have 



p K (x) = 1 - e~P tx 



14 
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with 
P= xp K (x)dp,(x) 



p K (x)dj2(x) 



-ptx^ 



To find the weight w(C\) of C\(H n ), we note that if f{x) : 
Z + ^R satisfies, using flM}, ± J\ /(&») = ~ Ei ^ = \F n \/n 
and thus [5|, Theorem 9.10] applies and yields 



dju(x). (4.13) 

:= x, then / : 
l = //dM, 



■;? 



-w(Ci 



n ^^ 



xp K (a;)dju(a;) = p. 



Combining f|4.13|) and (|4.14p . we thus find that 

\C 1 (G(n,t;F n ))\=w(C 1 (H n )) = pn + o p (n), 

where p solves the equation (|4.13p . which also can be written 

p = E(l-e- ptz ) = l-Ee~ ptz . 

Applying Lemma |4. II to Y := tZ, we see that when t > l/s2 
there is a unique p > satisfying (|4.16|) . 



(4.14) 

(4.15) 

(4.16) 

1/EZ, 



Further, in (ii) we may apply Lemma |4. II also to Y := tZ n ; thus there 
indeed is a unique such p n . Moreover, by Lemma l4njTiv) p n — > p. Hence, 
(I4TT51 yields 

\Ci(G(n,t;F n ))\ = p n n + o p (n), 



which proves (ii) when t > l/s2- 



We have shown the conclusions in (i) and (ii) when t < l/s2 and t > 



1/^2, respectively. However, the statements use instead the slightly different 
conditions t < l/s2(F n ) and t > l/s2(F n ). For (i) this is no problem: if 



t < l/s2(F n ) for infinitely many n, then t < l/s2 since we have assumed 
S2(F„) -> s 2 . 



To complete the proof of (ii) however, we have to consider also the case 
1/S2 > t > l/s 2 (F n 
ts 2 (F n 



). If this holds (for a subsequence), then M(tZ n ) 
< S2{F n )/s2 —}■ 1, and thus p n 
t < 1/^2, (|4.12p applies and shows that 



\Ci(G(n,t;F n ))\ = o p (n) = p n n + o p (n) 



by Lemma HJl^iv) Since 

(4.17) 



so 



tS 2 



(ii) holds in this case too. This completes the proof of (i) and (ii) 

Ve have, by 
t 3 S4. Hence 

25 n S2 



(iii) now follows easily from Lemma 14.11 We have, by (12. 5p . ~E(tZ n ) 
l + 5 n s 2 ,E(tZ n ) 2 



E(tZ n ) - 1 



E(tZ n ) 



t 2 s 3 and E(tZ n y 

$nS2 _ 5 n S2 

t 2 S 3 (1 + 5 n S 2 ) 2 S 3 



„3 



S3 



(4.18) 



so the lower bound follows by (ii) and Lemma 14. lltii 



For the upper bound we have by (I4.18P 

E(tZ n 



E(tZ n 



1 s s 



S3 
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and similarly 

(K(tZ n ) - l)E(tZ n ) 3 _ 5 n s 2 t 3 s 4 _ 5 n s 2 2 s 4 5 n sls 4 



(E(tZ„) 2 ) 2 t4 s 2 (l + <5 nS2 ) s 2 s | 



and the upper bound follows by Lemma 14. l|(iir 



For (iv) we note that if lim inf n 5 n > 0, we can by ignoring some small n 
assume that inf n 5 n > 0, and then the difference between the left-hand side 
and right-hand side in (J4.18P is bounded below (since 1 < s 2 < S3 = 0(1)); 
hence we can add some small r/ > to the right hand side of (I4.18P such 
that the inequality still holds for large n. Consequently, 

s 3 
d(G(n, t; F n ))/n > 5 n ^{\ - 25 n s 2 ) - rj + o p (l), 

s 3 

which implies that w.h.p. 



Ci(G(n,t; F n ))/n > 5 n ^(l - 25 n s 2 ) 



S3 
The upper bound follows in the same way. □ 

5. Proof of Theorems I3.2H3.5I 

Proof of Theorem \3.2l Define the functions 

f(t) := l/s 2 (t), 

g(t) := s 3 (t)/4(t) = f(t)s 3 (t), 
h(t) := s 4 (t)/4(t) = f 4 (t)s A (t). 

The differential equations (|3.2p - (|3.4p then translate into, after simple cal- 
culations including some cancellations, 

f(t) = -x\(t)f(t)-(l-xl(t)), (5.1) 

g'(t) = 3xl(t)f 3 (t) - 3x 2 (t)/(t) ff (t), (5.2) 

ti(t) = 7xl(t)f 4 (t) + 3(1 - xl(t))g\t)f-\t) - Ax\{t)f{t)h(t). (5.3) 

Consider first (|5.ip . The right hand side is locally Lipschitz in t and /, 
and thus there exists a unique solution with /(0) = 1 in some maximal 
interval [0,tf) with tf < 00; if tf < 00 (which actually is the case, although 
we do not need this), |/(t)| — > 00 as t f- tf. Since < x\(€) < 1 for all 
t > 0, and further x\(t) is decreasing, f'{t) < —(1 — x\{t)) < — cq, for some 
Co > and all t > 0.1, say. Hence, f{t) decreases and will hit at some 
finite time t c < tf. This means that s 2 (t) = l//(t) — > 00 as t /* t c , so (|3.2p 
has a (unique) solution in [0, i c ) but not further. 

We have f(t c ) = and thus, by flEE}, f'(t c ) = -(1 - x 2 (t c )) < 0. Conse- 
quently, defining p := (1 — S 2 (t c )) _1 > 0, 

f(t) = p- 1 (t c - t) ((1 + 0(t c -t)), t<t c , 
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and thus 

Mt) = ^-^{(i + o(t c -t)), t<t c , 

as asserted. 

Next, treating x~i(t) and / as known functions, (|5.2p is a linear differential 
equation in g. An integrating factor is 

G(t) :=3 / xl(u)f(u)du, (5.4) 

Jo 

and then the unique solution in [0,tf) is given by 

g (t) = e~ G{t) + 3e- G(t) [ e G[u) xl{u)f{u) du. (5.5) 



JO 

Hence f)3.3j) has the unique solution 3(t)^(t), i G [0, t c )> with g(i) given by 
(1531) . Note that g(t) > for i < t c . 

Let /3 := g(t c ) > 0. By (|5T2"j) . g'(t c ) = 0, and thus, for t < t c , g(t) = 
(3 + 0(t c - t) 2 , and 

s 3 (t)=(3sUt)(l + 0(t c -t) 2 ) = ^-^((l + 0(t c -t)), 

Finally we consider (|5.3|) . Here the right-hand side is singular at t c because 
of the factor f~ 2 (t) in the second term, so we modify h and consider 

hi(t) := h(t) - 3g 2 (t)s 2 (t) = h{t) - 3g 2 (t)r\t), 

which satisfies the differential equation 

h[(t) = 7x 2 (t)f\t) - I8x 2 (t)g(t)f(t) + I5x 2 (t)g 2 (t) - Ax 2 {t)f{t)h(t) 

= 7x 2 (t)f\t) - 18x 2 (t)g(t)f(t) + 3x 2 (t)g 2 (t) - 4x?(t)/(t)fci(t). 

Again, this is a linear differential equation, with a unique solution in [0, tf). 
We leave the explicit form to the reader, since we need only that h\(t) = O(l) 
for t < t c , which yields that for t G [0,t c ), 

s A (t) = h(t)4(t) = 3g 2 (t)s 5 2 (t) + h^sKt) 

= 3/3 2 4(t) + 0(4(t)). □ 

Proof of Theorem \3.3[ For k = 2, this is, as said above, proved in [14J . The- 
orems 1.1 and 4.3]. We prove the extension by the same method (with 
somewhat different notation). 

Let, for a vertex v G G, C(v) be the component of G containing the vertex 
v, and C(v) := \C(v)\. 

For a given graph G, let G + be the random graph obtained by adding 
one random edge by the Bohman-Frieze rule; we assume that the edge was 
chosen from the pair e\ = {v\,w\} and e 2 = {^2,^2}- If the added edge is 
{v,w} (which thus is either {v\,wi} or {^2,^2}), and further C(v) 7^ C(w), 
then, by (l2TT|) . 

S k (G + ) - S k (G) = (C(v) + C{w)) k - C{v) k - C{w)\ (5.6) 
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while S k {G + ) - S k (G) = if C(v) = C(w). We define 

A* k = A* k (G; v, w) := (C(v) + C(w)) k - C{v) k - C(w) k . (5.7) 
Hence, 

E(S k (G+) - S k (G) - A* k ) = -E(A* k l[C(v) = C(w)]) 

= -E((2 fc - 2)C(v) k l[C(v) = C(w)]) 
and thus 

| E{S k (G + ) - S k (G) -A* k )\< 2 k E(C(v) k l[C(v) = C(w)]) 

< 2 k E{C(v 1 ) k l[C(v 1 ) =C(w x )]) +2 k E{C(v 2 ) k l[C(v 2 ) =C(w 2 )]) 

= t EC( Vl ) k+l < 2 k+1 ^> . 

n n 

In particular, if C\{G) = O(logn), then 

| E{S k {G + ) - S k (G) -A* k )\ = o(^^) = o(l). (5.8) 

Expanding (|5.7p . we have 

A^ = 2C{v)C(w), (5.9) 

A* 3 = 3C{v) 2 C(w) + 3C(v)C(w) 2 , (5.10) 

A^ = AC{v) 3 C(w) + 6C(v) 2 C{w) 2 + 4C(v)C{wf. (5.11) 

The Bohman-Frieze rule is to take {v,w} = {vi,w\} if C{v\) = C(w{) = 
1. The probability of this is x\(G) 2 , and in this case A* k = 2 k — 2. 

The opposite case {v,w} = {^2,^2}, which we denote by £2, has proba- 
bility 1 — x\{G) 2 . Conditioning on this case places us basically in the well- 
studied Erdos-Renyi regime. That is, v and w are uniform and independent, 
and thus for any k and £, 

E( C (v) k c( W y 1 £2) = ± e C (v) k cH e = ± E c " +l E c r 

v,w i j 

= s k +i(G)si +1 (G). 
Hence, (pT9l ~ ([5TLT) yield 

E A; = 2x\{G) + (1 - Xl (G) 2 ) ■ 2s 2 (G) 2 , 
E A* 3 = Qx\{G) + (1 - xi(G) 2 ) • 6s 2 (G)s 3 (G), 
EAl = Ux 2 (G) + (1 - x^Gf) ■ (8s 2 (G)s 4 (G) + 6s 3 (G) 2 ). 
By (|5.8p . we thus have, for k = 2, 3,4 and provided C\{G) = O(logn), 

E(S k (G + )-S k (G))=EA* k + 0{log k+1 n/n) 

= 2f k (x 1 (G),S2(G),s 3 (G),s±(G))+0(log k+1 n/n), (5.12) 
with 

/ 2 (xi, 82, s 3 , s 4 ) := x\ + (1 - x^s 2 ,, 
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f3(xi,s 2 , s 3 , Si) := 3x\ + 3(1 - xl)s 2 s 3 , 

/ 4 (xi,S2,s 3 ,s 4 ) := lx\ + (1 -xf)(4s 2 s 4 + 3s|). 

Similarly, as shown in [141 ] . 

E(n 1 (G + )-n 1 (G))=2f 1 (x 1 (G),s 2 (G),s 3 (G),s A (G))+0(l/n), (5.13) 

where (the variables S2,S3,s 4 are redundant here) 

fi(xi,s 2 , s 3 , s 4 ) := -x? - (1 - xf)xi. 

Consider the vector-valued random process 

Xi := (x 1 (BF i ),s 2 (BF l ),s 3 (BF i ),s 4 (BF l )), 

and let F = cr(Ao, . . . , Xi) be the cr-field describing the history up to time 
i. Further, let $ := (f 1} f 2 , / 3 , / 4 ) : R 4 -> M 4 . Using this notation, ([5^2]) - 
(I5TT3]) yield 

E(n(X i+1 - Xi) | 7=i) - 2$(Xi) = 0(log 5 n/n), (5.14) 

uniformly in i < in/2, provided Ci(BFj) = O(logn). 

By [14, Theorem 1.1], there exists a constant c' (depending on t) such that 
w.h.p. Ci(BFj) < c'logn for all i < in/2. As in [1J], we avoid the problem 
when Ci(BFi) > c'logn by defining X * = X = (1,1,1,1), X* +1 = X m 
when Ci(BFj) < c'logn and X* +1 = X* + |$(X*) otherwise. Then w.h.p. 
X* = Xi for all i < in/2, so we can just as well consider X*. We have, by 
(|5.14p but now without side condition, for all i < in/2, 

E(n(X* +1 - X*) | Fi) = 2$(A7) + 0(log 5 n/n) 

and also, for some c", from (|5.6p and \n\(G + ) — m(G)\ < 2, 

\X* +1 - X* | <c" log 4 n/n. 

The differential equation method i n the form of ISpencer and Wormaldl 14 . 



Theorem 4.1], which is taken from IWormaldl la . Theorem 5.1], now applies 



(with Y(i) = nX*) and the result follows; note that the differential equations 
(|3.ip - (l3.4p can be written <p'(t) = $((^(i)) with (p = (x\, s 2 , s 3 , s 4 ), where 
further ip(0) = (1,1,1,1) = X = X*. D 

Proof of Theorem \3.5[ We may assume that 5 is small, since the result is 
trivial for 5 > 5o > if we choose K large enough. In particular, we assume 
5 < 1. 

Let e := <5 2 / 3 > 5. We stop the process at i c — e, and let F := BF(i c — 
e). We then let the process evolve to i c + 5 by adding (e + 5)n/2 further 
edges according to the Bohman-Frieze rule. Actually, for convenience, we 
add instead a random number of edges with a Poisson distribution Po((e + 
5)n/2) ; this will not affect our asymptotic results (by the same standard 
argument as for comparing the different models in Section [2]) . We denote 
the resulting graph by BF(i c + 5). 
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By Theorems 13.31 and 13, 2\ for k = 2,3,4, and with a*; as in Theorem 13.21 
s k (F) = s k (t c -e) + o p (l) = -^3 (1 + 0(e)) + o p (l). 



Since |o p (l)| < e w.h.p., we thus have w.h.p. 

^(F) = J^(l + 0(e)). (5.15) 

(This means that there exists a constant c, not depending on e or n, such 
that (|5.15|) holds with the error term 0(e) £ [—ce,ce] w.h.p.) Similarly, 

xi(F) -?-> xi(t c - e) = xi(t c ) + 0(e), so w.h.p. x x (F) = x\{t c ) + 0(e). 

We fix F (i.e., we condition on F) and assume that (J5.15P holds together 
with x\(F) = x\(t c ) + 0(e) (for some fixed implicit constant c in the 0(e); 
we have just shown that this holds w.h.p. provided c is chosen large enough). 

We cannot directly apply Theorem 12.11 since the graph evolves by the 
Bohman-Frieze evolution and not by the Erdos-Renyi evolution. Neverthe- 
less, we can approximate and find upper and lower bounds of the graphs 
where we can apply Theorem 12.11 the idea is that we consider the Erdos- 
Renyi edges separately as an Erdos-Renyi evolution. 

For a lower bound, let V\ be the set of isolated vertices in F and consider 
only the pairs of edges e± = {vi,w\}, e 2 = {v2, ^2} where v\ £ V\ or w\ ^V±. 
Since the graphs BF^ in the continued process contain F, the vertices v\ and 
w\ are not both isolated in the current BF^, and thus e2 = («2j ^2) is added, 
and these are independent Erdos-Renyi edges, i.e., uniformly chosen. The 
number of such Erdos-Renyi edges is Po((l — x\(F) 2 )e + S)n/2) , since each 
time we add an edge, the probability of it being of this type is 1 — (|Vi|/n) 2 = 
1 — x\(F) 2 . (Note that we ignore some Erdos-Renyi edges in order to avoid 
unpleasant dependencies.) 



Call the resulting graph H C BF(t c + J). Then Theorem [271 If iv) applies 
to H~ , with 

t = (1 - x 1 (Ff)(e + 5) = (1 - xi(t c ) 2 + 0(e))(e + 5) 
and, recalling (|5.15p and a = (l — xf(t c )) , 

5 n = t-l/s 2 (F) = (l-x 1 (t c ) 2 )(e+S)-a- l e+0(e 2 ) = (l-x 1 (t c ) 2 )5+0(e 2 ), 

(5.16) 
which yields w.h.p., using again (|5.15p . 

Ci(BF(* c + £)) ^ C X (H-) ^ s 2 (Ff H 

> > M n 7=r-(l - 26 n S2(F)\ 

n n S 3\F) 

= 2 ((1 - Xl (t c ) 2 )8 + 0(e 2 )) ^, (l + 0(e) + O (^ 

= ^5(1 + 0(e) + 0(5/e)+0(e 2 /5)) 

= ip 6 ( 1+ ° {6l/3) ) =76+ 0(54/3) ' ( 5 - 17 ) 
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with our choice e = <5 2 ' 3 (which is optimal in this estimate). 

For an upper bound, note that w.h.p. at most (e + 8)n < len edges are 
added to F, so at most Aen vertices are hit, and thus during the process 
from F to BF(t c + 5), 

x\ > xi(F) - 4e = xi (t c ) - 0{e). 

Hence we add w.h.p. at most 

(1 - (xi(t c ) - 0{e)f){e + 5)n/2 = (l - Xl (t c ) 2 + 0(e)) (e + <5)n/2 

Erdos-Renyi edges. We also add a number of non-Erdos-Renyi edges, all 
joining two isolated vertices (or being loops). They may depend on the 
Erdos-Renyi edges already chosen, but we avoid this dependency by being 
generous and adding the edge e\ = (ui,itfi) in each round whenever both 
v\ and w\ are isolated in F and neither is an endpoint of an already added 
non-Erdos-Renyi edge. (We add e<i by the same Bohman-Frieze rule as 
before, so we may now sometimes add both e\ and e 2 -) 

Let c\ be a large constant and let H + be the graph obtained from F 
by adding 2en (to be on the safe side) non-Erdos-Renyi edges in this way, 
together with {l — x\ (t c ) 2 + c\e\ (e + 5)n/2 Erdos-Renyi edges, independent 
of each other and of the non-Erdos-Renyi edges. We conclude that, if c\ is 
chosen large enough, we may couple H + with the Bohman-Frieze process 
such that w.h.p. BF(t c + 8) C H+. 

Since the two types of edges are added independently, we may further add 
all non-Erdos-Renyi edges first. Let F\ be F together with all non-Erdos- 
Renyi edges. There are 2sn such edges, and each joins two isolated vertices 
and changes Sk by 2 k — 2 (or by if the edge is a loop). Hence, for every 
k < 4, by fl5J5]) . 

SfcCFi) = s k (F) + O(e) = ^ (1 + 0(e)) . (5.18) 



Since H + is obtained by adding the Erdos-Renyi edges to F\, Theorem l2.ll 
applies with 

t= (l- Xl (t c ) 2 + Cl e)(e + 5) 
and 

5 n = t-l/s 2 (F 1 ) = (l-xi(t c ) 2 )(e+5)-a- 1 e+0(e 2 ) = (l-xi(£ c ) 2 )£+0(e 2 ), 

(5.19) 
the same estimate as was obtained in (|5.16p . We use the upper bound in 



Theorem Ejiv)| By (^T8l) and (^T9l) . 



^^"fl^^^T)^©^^' 
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Hence Theorem [27l Ifiv) applies (for small 5) and yields, w.h.p. 

S2(F 1 ) 3 



n^dCBFfo + 5)) < n- 1 C l {H + ) < 25 n -^-{l + O^)) 

a 3 e~ 

f3a 3 e 



2(1 - Xl (t c ) 2 )d-^-^ (1 + 0(e 2 /5 + e + 5 1 / 3 )) 



= 2(1 ^ l(tc) ^ (1 + O^ 1 / 3 )) = 7 J + 0(5 4 / 3 ). 
This and the corresponding lower bound (I5.17P yield the result. D 
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